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Abstract 

We develop a class of soliton solution of linear Schrodinger equation with- 
out external potential. The quantum probability density generates its own 
boundary inside which there is internal vibration whose wave number is de- 
termined by the velocity of the particle as firstly conjectured by de Broglie. 
Assuming resonance of the internal vibration will lead to quantization of par- 
ticle's momentum in term of wave number of the envelope quantum proba- 
bility density. We further show that the linearity of the Schrodinger equation 
allows us to have non-interacting many solitons solution through superposi- 
tion, each describing a particle with equal mass. 
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1. Introduction: de Broglie's double solution 

In his attempt to explain the dual nature of matter as both wave and 
particle, de Broglie proposed what is called as the theory of double solution 
He was searching for a nonlinear type of wave equation which possesses 
a singular, spatially localized, nondispersive solution in region where the 
amplitude is sufficiently high. This spatial part of the global solution is then 
supposed to play the role of a particle. On the other hand, in region where 
the amplitude is weak, the solution is assumed to behave like a linear wave 
obeying the law of superposition. He thus envisioned a particle which is 
immersed in linear wave field and is guided by the latter. He argued that in 
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order that each point of the singular solution to move non-dispersively, thus 
particle-like, the velocity of the singular solution must be equal to the velocity 
of the non-singular linear solution. This is de Broglie's guiding principle |l| 
which eventually led him to the famous relation between the momentum of 
a particle p and the associated wave length A of the guiding wave: 

P — hk, (1) 

where k = 2tt/X is the corresponding wave number. 

The prevailing quantum theory is however founded on the Schrodinger 
equation which is linear with respect to its ingredient, namely the wave 
function. In this theory, a single free particle is usually represented by a 
plane wave. Yet, the successes of this representation seems to rely on the 
adoption of Eq. ([1]) where k is now the wave number of the plane wave, 
and the relation between the energy of the particle and the frequency of the 
plane wave, to: E = hw One can even argue that those two relations 
are the most important principles of quantum theory 0, 0] . Even the linear 
Schrodinger equation can be derived from those relations. 

It is then apparently the linearity of the Schrodinger equation which dis- 
courages people to keep searching for particle-like solution. The ambition to 
search for a particle-like solution is further overshadowed by the tremendous 
successes (for all pragmatical purposes 0) of the probability interpretation 
of wave function invented by Born [g| which is heavily based on the law of 
superposition. This law is possible only if the underlying equation is lin- 
ear. This eventually led to the axiomatic development of quantum theory by 
Dirac and von Neumann 0, 0] . 

Despite of these facts, fueled by the dissatisfaction on the foundation of 
the axiomatic approach to quantum theory 0, 0], and a simple fact from, 
say the single slit experiment which undeniably suggests that there is indeed 
a localized quantity which hits the screen thus giving a point of scintillation 
[1, Gil > some people are still searching for a quantum theory which assumes 
a localized particle-like solution. Interestingly, this adventure is again led 
by de Broglie after abandoning his own idea for many decades. Much effort 
again is made to search for such a solution from a nonlinear wave equation 
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In this paper, we shall show that even a linear Schrodinger equation pos- 
sesses a class of localized-nondispersive solution which interestingly shares 
the properties long conjectured by de Broglie. This will be done by enter- 
taining the hydrodynamics picture of the Schrodinger equation developed by 



2 



Madelung 15|. Even more, we show that the soliton wave function is the 
limiting casejrf a class of the most probable wave functions given its quan- 
We then show that the linearity of the Schrodinger equation 
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turn energy 

allows one to develop many solitons solution through superposition. 



2. Most probable quantum probability density 

Let us consider the following linear Schrodinger equation which is sup- 
posed to describe a single free particle with mass m 
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iW^{q ] t) = -^S^{q ] t). (2) 
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Here q is space and t denotes time. For our discussion through out this 
paper it is sufficient to consider spatially one dimensional space. i>(q;t) 
is a complex-valued wave function supposed to determine the state of the 
particle. Let us project the above dynamics into real space. To do this, let 
us make the following transformation if)(q; t) = R(q; t) exp(iS(q; t)/h), where 
R(q] t) and S(q; t) are real-valued functions. Putting this into Eq. (jSJ) and 
separating into the real and imaginary parts, one obtains 

dv 

m— = -d q U, d t p + d q {yp) = 0. (3) 

Here, p(q) = |V ; ( ( ?)| 2 = R 2 (q) is quantum probability density, U(q) is deter- 
mined by the quantum amplitude R(q) as 

h 2 d 2 R 

™ ^ -srJr (4) 

and v(q) is a velocity vector defined by the quantum phase S(q) as 

v{q) = d q S/m. (5) 

Equation (EJ) to (EJ) is the so-called Madelung fluid picture of the Schrodinger 
equation [3]. Due to the formal similarity with Euler equation, the term 
on the right hand side of the left equation in Eq. F = —d q U, is called 
as quantum force, thus accordingly, U(q) is called as quantum potential. 
Putting in Madelung fluid form, it thus became explicit that the original 
Schrodinger equation possesses a hidden self-referential property. Namely, 
U(q) is determined by p(q) and in turn U(q) will dictate the way p(q) evolves 
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with time through Eq. ([3]) and so on and so forth. One thus may expect to 
observe self-organized physically interesting phenomena. 

Next, let us consider a class of quantum probability densities which max- 
imizes Shannon information entropy: H[p] = — J dq p(q) In p(q) [T3], given 
its average quantum potential U = J dq U(q)p(q). This is the so-called max- 



imum entropy principle [18j . It has been argued as the only way to infer 



from an incomplete information which does not lead to logical inconsistency 



191 ]. In our present case, the limited information that we have in hand is the 
average quantum potential. We shall show later that the average quantum 
potential is a physically relevant information. Hence, the maximum entropy 
principle will give us the most probable quantum probability density with 
average quantum potential U . This inference problem can be directly solved 
to give 2o| : 

p(q) = ^ex V (-U(q)/T), (6) 

where T is a constant determined by U, and Z(T) is a normalization factor. 
Combining with Eq. (TjJ, Eq. (jHJ) comprises a differential equation for p(q) 
or U(q), subjected to the condition that p(q) must be normalized. In term 



of U(q), one has [16 



1 AmT 

dlu = ^{d q uf + ^u. (7) 

Figure [T^i shows the solution of Eq. (J7|) with the boundary conditions: 
£7(0) = 1, d q U(0) = 0, for several small values of positive T. The reason 
for choosing small values of T will be clear as we proceed. All numerical 
results in this paper are obtained by setting h = m = 1. One can see that 
the quantum probability density is being trapped by its own self-generated 
quantum potential. Moreover, one can also see that there are points q = 
±L m , where the quantum potential is blowing-up, namely f/(±L,) = oo. 
This is a familiar phenomena in nonlinear differential equation [2l| and in 
particular, for the case we are considering, it can be proven as follows. Let 
us define a new variable u(q) = d q U . The nonlinear differential equation of 
Eq. (jZ|) then transforms into 

1 , AttiT . . 

B * u= vf* + -iF u - (8) 

The boundary condition translates into u(0) = d q U(0) = 0. Further, let us 
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now consider the following nonlinear differential equation 

^ = ^ + (9) 

where X = £7(0); with u(0) = 0. Since U(q) > £7(0) = X, then it is obvious 
that > |w(<?)|- 




Figure 1: (a) The quantum probability density (solid line) and its corresponding quantum 
potential (dashed line) for several small values of T obtained by solving equation jJJ . We 
also plot the analytical solution for po(q) at T = 0, assuming that Uo(q) takes the form of 
a box with infinite wall. See text for detail, (b) u(q) and u(q). See text for detail. 

On the other hand, one can solve the latter nonlinear differential equation 
of Eq. Q analytically to have 

2T f 

u(q) — atan(feg), a = —y2mX, b = —v2mX. (10) 

a h 

It is then clear that at q — ±q = ±7r/ (26), u is blowing-up, namely u(±q) = 
±oo. Recalling the fact that \u(q)\ > (^(q 1 )!, then u(q) is also blowing-up at 
points q = ±L m , u(±L m ) = ±oo, where L m < q. See Fig. [lb. Hence, one 
can conclude that U(q) is also blowing-up at q — ±L m , U(±L m ) = oo. It is 
then safe to say that the self-trapped quantum probability density possesses 
only a finite range of spatial support M. = [—L m , L m ]. 

In Fig. [21 we plot the variation of the half length of the support against 
T. One observes that as T increases toward infinity, L m (T) decreases to- 
ward zero. This shows that the corresponding quantum probability density 
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Figure 2: The half length of the support of quantum probability density plotted against 
the parameter T . 

collapses toward a delta function. On the contrary, as one decreases T to- 
ward zero, we see from Fig. [2] that L m is approaching to a finite value, 
lim r ^ L m = Lq. This hints us that the quantum probability density and 
thus its corresponding quantum potential are also converging toward some 
functions 

KmU(q;T) = Uo(q), lim p(q; T) = p (q). (11) 

To discuss this latter situation in more detail, let us go back to Fig. UK 
One can see that as T decreases, the quantum potential inside the support 
is getting flatterer before becoming infinite at the boundary points, q = 
±L m (T). One can thus guess that at vanishing T, T = 0, the quantum 
potential is converging toward a box of length 2L with perfectly flat bottom 
and an infinite wall at the boundary points, q = ±L . Let us show that 
our guess is correct. To do this, let us proceed to calculate the profile of 
the quantum probability density for vanishing T, po(q), assuming that U (q) 
takes a box potential of length 2Lq. Namely Uo(q) = U c , where U c is constant 
for — L < q < L , and Uq(±L ) = oo. Since Uo(q) is flat inside the support, 
then one has Uq = J^l dqp (q)Uo(q) = U c . Hence, from the definition of 
quantum potential given in Eq. (TJJ, inside the support of the quantum 
probability density, one has 

- ^-d 2 q R = U C R = U Ro, (12) 
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where Ro(q) = pQ 2 (q) is the quantum amplitude at T = 0. The above differ- 
ential equation must be subjected to the boundary condition: i? (±L ) = 0. 
Solving Eq. f|T21) one obtains 



Ro(q) = A cos(k q) 



(13) 



where A is a normalization constant and k is related to the average quantum 
potential as 



The boundary condition imposes /co-^o — 7r/2. In Fig. [Th, we plot the 
above obtained quantum probability density, po(q). One can see that as T is 
decreasing toward zero, p(q; T) obtained by solving the differential equation 
of Eq. ((7j) is indeed converging toward po(q) given in Eq. ([TBI . This confirms 
our initial guess that at T = the quantum potential takes a form of box 
with infinite wall at q — ±L . 

Let us remark that all the above results suggest that the wave function can 
generate its own boundary at which it is vanishing. This is a very important 
fact since then one can talk about a self-sustained geometrical property of 
the wave function. We shall show in the next section that this fact will help 
us in developing our notion of particle from the Schrodinger wave function. 
One can of course trace back the emergence of this self-organized structure 
from the self-referential property of the Madelung fluid dynamics. 

3. Soliton wave packet 

Let us now choose a pair of fields {po(q) , v (q)} as the initial state of the 
Madelung fluid. Here, po(q) = Rl(q) is given in Eq. (fTBj) and v (q) = v c 
is a uniform velocity field with non-vanishing value only inside the spatial 
support M. = [—Lq,Lq]. Then, since at t — the quantum potential is flat 
inside the support, the quantum force is initially vanishing, d q U = such 
that dv/dt = 0. Hence, the velocity field at infinitesimal lapse of time t = At 
will stay unchanged and keeps uniform. This in turn will shift the quantum 
probability density in space by Aq = v c At while keeps its profile unchanged. 
Correspondingly, it will shift the support of the quantum probability density 
by the same amount. The same thing occurs for the next infinitesimal lapse of 
time and so on and so forth. One thus concludes that the quantum probability 
density is moving with a uniform velocity field v c , keeping its initial form 




(14) 
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remained unchanged. Hence, at time t, the quantum probability density is 
given by 



p(q; t) = p (q - v c t; 0) = Al cos 2 (k q - u t), (15) 

where Uo = koV c and q G M t = [v c t — L , v c t + L ]. 

To write the explicit form of the complex- valued self-trapped and uniformly- 
moving yet stationary wave function i/j(q; t) of our single particle, let us first 
calculate its quantum mechanical energy. Since quantum mechanical energy 
is conserved then it is sufficient to use the wave function at t — 0. Putting 
the wave function in polar form, ijj(q) = Ro(q) exp(iS(q)/H), one has 

-^Rod 9 Rod q S-^^ q S). (16) 

The first term on the right hand side is equal to the average quantum poten- 
tial, Uq — J dql/oPo = U c = h 2 k\j{2m). Next, defining kinetic energy density 
as K (q) = (m/2)t>Q(g), then the second term is equal to the kinetic energy of 
the Madelung fluid K = J dqK (q)p (q) = (m/2)v 2 c . Further, for a uniform 
velocity field, the last term is vanishing, (l/m)d 2 S = d q Vo = 0. For the same 
reason, since Ro(q) is an even function and d q R (q) is an odd function then 
the third term is also vanishing. 

Hence, in total, the quantum mechanical energy of the self-trapped wave 
function moving with a uniform velocity field v c can be decomposed as 

<S>=». + *. = ^ + ==£ (17) 

Moreover, using similar argument as above, the average quantum momentum 
can also be calculated to give 

(p) = / ° dqil)*(q)(-ihd q )ij)(q)=mv c . (18) 

J-Lo 

The average kinetic energy of the Madelung fluid and average quantum mo- 
mentum are thus related as K = (p) 2 /(2m). This observation leads us to 
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conclude that Uq must be interpreted as essentially an internal energy of the 
single particle. Namely it is the energy when the particle is not moving. 
A similar notion is also developed in special relativity theory, dubbed as 
rest mass energy. One can also check that the internal energy is the energy 
which is missed to be taken into account if one uses the simple plane wave 
to represent a free moving particle. 

Bearing in mind the above derived facts, let us write the moving- stationary 
pair of fields {p(q), v(q)} in complex- valued form, ip(q) = R(q; t) exp(iS(q; t)/K) 
To do this, one has to calculate the quantum phase S(q; t) by integrating 
v(q) = {l/m)d q S = v c to give S(q;t) = mv c q + £(t), where £(t) depends only 
on time t. The wave function we are searching thus takes the following form: 

rfsiQit) = A cos (ko(q — i> c t))exp (i(mv c q + ^(t))/h\. Finally, putting this 

back into the Schrodinger equation of Eq. fl2j), keeping in mind u = k v 
and Eq. (IT7j) . one easily obtains d t £ = —(E), which can be integrated to give 
£(t) = — (E) t modulo to some constant. We are thus led to the following 
form of wave function: 



where q G M. t = [(p)t/m — L Q , (p)t/m + L ]. The solitonic nature of the 
above wave function is obvious. 

4. Compton wave length, phase resonance, and superposition 

Let us proceed to discuss the stationary-moving soliton solution of the free 
particle Schodinger equation that we just developed in the previous section. 
First, the wave function possesses only a finite range of support Aq = 2Lq. 
Using the fact that koLo = tt /2 which came from the boundary condition 
for the quantum amplitude Rq(L ) = 0, and the expression for k in term 
of quantum potential U c = U of Eq. (JHJ), one is led to the following rela- 
tion between the width of the support of the soliton wave function and the 
quantum potential: 



On the other hand, we have shown in the previous section that Uq is but 
the internal energy, namely the energy of the particle when it is not moving. 
It is thus reasonable to assume that this energy is equivalent to the rest mass 




(19) 




(20) 
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energy predicted by the Einstein's theory of special relativity. Namely, one 
assumes that 

U = mc 2 , (21) 
where c is the velocity of light. Putting this into Eq. (I20I) . one finally obtains 

A ^7i-^7f Ac ' (22) 

where Ac = h/ (mc) is but the Compton wave length. One can thus conclude 
that the active region of the soliton wave function is proportional to the 
Compton wave length of the particle. In particular, this active region is 
narrower for particle with larger mass. 

Next, in Fig. [3] we plot the snapshot of spatial profile of the real part 
of soliton wave function given in Eq. ffl9|) . ip s (q;t), at times: t < t 1 < t 2 . 
One observes an envelope given by Ro{q — v c t) which is moving uniformly 
with velocity v c . One can also see from Eq. ( fl9i) . that there is internal spa- 
tial vibration whose wave number ks is determined by the average quantum 
momentum 

^<|U- (23) 

\b = 2-K/ks = 27rh/(mv c ) is then but the famous de Broglie wave length. 
This is in fact what is expected by de Broglie in his conjecture of double 
solution. Yet in contrast to the theory of double solution which envisions 
a soliton solution emerging from a nonlinear wave equation, here we have 
shown a soliton wave function as a solution of linear Schrodinger equation. 

Before proceeding to discuss further consequence of the linearity of the 
Schrodinger equation, let us evaluate what will happen if the internal vibra- 
tion of de Broglie's wave length resonates inside the self-generated quantum 
potential. Namely one imposes that the length of the support of the quan- 
tum probability density, Aq = 2L , is equal to the integer multiple of the de 
Broglie's wave length X B . Hence, k B is equal to the integer multiple of k , 
ks = nko, where n = 0, ±1, ±2, .... In this case, the velocity of the particle 
and the wave number of the envelope quantum amplitude ko are related as 

(p) = mv c = nhk , n = ±1, ±2, .... (24) 

Namely, in this case the momentum of the particle can only take discrete 
possible values with the spacing given by po = hk . In the next section, 
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Figure 3: The snapshot of the real part of soliton wave function, ip s (q;t), at t Q < t\ < t%. 



we shall discuss a possible physical situation which requires the above phase 
resonance thus induces quantization. 

Finally, let us discuss important consequence of the linearity of the Schr6ding< 
equation. Equation (fl9l) tells us that given the energy (E) defined as in Eq. 
( TTTjl and the momentum (p) of the single particle of mass m, one can associate 
a single soliton occupying a localized space moving with velocity v = (p)/m. 
Since the Schrddinger equation is linear, then any superposition of many soli- 
ton solutions of the type given in Eq. / fTPj) will also satisfy the Schrddinger 
equation. For example, let us consider the following superposition of two 
solitons solution 



Each soliton term on the right hand side moves with uniform velocity 
i — 1, 2, and possesses a finite support given by M. i = fat — L^v^t + Li], 
where Li = 7r/(2ki). The total support is then given by M. ms = M.\ © 
One can then show that the above wave function satisfies the Schrddinger 




(25) 



equation of Eq. fl2]) if Ei is given as follows: 

E ' = l^ + 2 mV >- i = 1 ' 2 ' (26) 

Hence, it is as if each term on the right hand side of Eq. ( 1251) describes a 
de Broglie's particle of mass m with internal energy Ui = h 2 kf/(2m), moving 
with velocity Vi thus possessing kinetic energy Ki = (l/2)mvf. Both of the 
particles are not interacting. Without loosing generality, let us consider the 
case when the sign of v\ and v 2 are opposite to each other. First, at t = 0, one 
observes spatial interference pattern. Hence, there is no way to distinguish 
one from the other. One thus must consider both as a single particle of 
mass m with energy (E). Yet, at sufficiently large time, the two solitons 
are spatially separated, thus no interference phenomena is seen. One can 
distinguish two spatially localized solitons corresponding to two particles. 
At this situation, it is easy to calculate the conserved quantum energy (E) 
to have 

(E) =E X + E 2 . (27) 

Namely, it is given as the summation of the energy of each particle, E^. The 
theory then predicts an interesting phenomena in which a particle of a given 
mass m can break into two particles with still equal mass m yet possessing 
different energy. Since the Schrodinger equation is time reversal, then by 
reversing the time, one concludes that two particles with equal mass m can 
merge into a particle of the same mass yet with energy given by the addition 
of the energy of the merging particles. 

One can of course extend the above two solitons solution to many solitons 
solution straightforwardly, each of which can be interpreted as describing a 
single particle of equal mass. Moreover, each moves independently from the 
other thus the total energy is given by the summation of the energy of all 
the particles. Hence, the Schodinger equation of Eq. (J2]) which previously 
is thought of as a theory of a single free particle, now is re-interpreted as a 
theory of non-interacting many particles with the same masses. 

5. Conclusion and Discussion 

First, we have shown that the linear Schrodinger equation for a single free 
particle possesses a class of soliton solution parameterized by the quantum 
energy and momentum. We showed that the soliton wave packet shares the 
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properties that a quantum particle should have as envisioned by de Broglie 
long time ago in the early days of quantum mechanics even before the birth 
of the Schrodinger equation. Namely, the soliton wave function possesses an 
internal vibration whose wave number ks is determined by the momentum 
of the particle (p) exactly described by the famous de Broglie's relation, 
(p) = hk B . 

We then showed that assuming resonance of the internal vibration inside 
the quantum potential, the possible values of the momentum of the particle 
is discretized as the integer multiple of the wave number of the envelope 
soliton. A natural question then arises: when this quantization happens? 
To answer this question, one can borrow conventional idea from the physics 
of classical wave. In this field, a resonance usually is related to interaction 
and transfer of energy. We therefore expect that the quantization in soliton 
wave function happens if it interacts with other physical object. We shall 
elaborate this idea as future work. 

Finally, we showed that the linearity of the Schrodinger equation allows 
us to construct many solitons solution through the superposition principle. 
We argue that each term of the superposition can be interpreted as describing 
a single particle with the same mass as the other, yet each might possesses 
different momentum and energy. Moreover, each particle is free to move 
independently from the other. Hence, we conclude that the Schrodinger 
equation of Eq. ([2]) which is conventionally thought of as a theory describing 
a single free particle, can now be regarded as a theory of many particles with 
equal mass. 

It is then interesting to further ask: can we have a similar theory which 
possesses many solitons solution so that each can be associated to a particle 
with different masses? To do this one can not start from a theory in which 
the mass is already fixed as in the Schrodinger equation of Eq. (TS]). Namely, 
one has to develop a theory in which the masses emerge as consequence of 
the theory. This theory should be based on a mass-less wave equation. We 
shall report this some where else. 
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